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Abstract 

The percolation transitions on hyperbolic lattices are investigated numerically using finite-size 
scaling methods. The existence of two distinct percolation thresholds is verified. At the lower 
threshold, an unbounded cluster appears and reaches from the middle to the boundary. This 
transition is of the same type and has the same finite-size scaling properties as the corresponding 
transition for the Cayley tree. At the upper threshold, on the other hand, a single unbounded 
cluster forms which overwhelms all the others and occupies a finite fraction of the volume as well 
as of the boundary connections. The finite-size scaling properties for this upper threshold are 
different from those of the Cayley tree and two of the critical exponents are obtained. The results 
suggest that the percolation transition for the hyperbolic lattices forms a universality class of its 
own. 
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I. INTRODUCTION 



Geometry is often crucial in physical phenomena, since dimensionality and topological 
defects determine the properties of phase transitions. Progress in complex networks is an 
example where geometrical features are paramount to the physical properties (for a review, 
see, e.g., Ref. This is also reflected in ongoing research interests in the nontrivial lattice 
structures such as fractal lattices and Apollonian networks j^]. Another important area 
is the hyperbolic geometry where even the familiar standard physical models like electronic 
spins and Brownian motions exhibit novel behaviors [4]. Such studies are interesting not 
only from a theoretical viewpoint but also from a potential real applicability related to the 
rapid development of fabrication of devices and structures in nanoscales 5|. In this work, 
we focus on the percolation problem and investigate the percolation phase transition for 
negatively curved hyperbolic lattices. Percolation on hyperbolic lattices has so far been of 
predominantly mathematical interest and several interesting mathematical results have been 
reported, including the existence of an intermediate phase with infinitely many unbounded 
clusters p, 7, sj]. More specifically, there exist in general two critical thresholds: unbounded 
clusters are formed at the first threshold, while unbounded clusters are merged to become one 
unique unbounded cluster at the second threshold. In case of the usual flat lattices, which 
have vanishing surface- volume ratios in infinite- volume limit, the two thresholds coincide, 
and the theory of such percolation transitions is well-developed [3]. On the other hand, 
the properties of the percolation transitions for hyperbolic lattices still remain to be further 
clarified. In particular, these lattices are not homogeneous due to the presence of a boundary 
and as a consequence, the critical properties may differ from the mean-field-type transitions 
which were discussed in Ref. 1(3]. In this paper, we investigate the characteristic features of 
the thresholds and the corresponding phases by using various statistical measures like the 
number of boundary points connected to the middle, the ratio of the first and second largest 
clusters, and the cluster size distribution. We use finite-size scaling methods to obtain the 
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critical properties, together with the Newman-Ziff algorithm 

This paper is organized as follows: In Sec. [TT1 we introduce alternative manifestations of 
the percolation thresholds and show that all of them coincide in an ordinary square lattice. 
In Sec. IHIt we explain the notion of a hyperbolic lattice and the appearance of two distinct 
percolation transitions. We start with the Bethe lattice and the standard mean-field results. 
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Next the Cayley tree is introduced which is the simplest model with two thresholds. This 
model is used as a benchmark when discussing the characteristic features of percolation in 
hyperbolic lattices. The efficiency of the finite-size scaling methods is also illustrated for 
this simpler model. Then we present the results for the general hyperbolic structures which 
contain loops. We summarize our results in Sec. IIVI 



II. VARIOUS MANIFESTATIONS OF PERCOLATION 

There are various possible manifestations of percolation. We here describe the critical 
thresholds corresponding to them as well as their relationship. The percolation threshold 
is usually defined as the occupation probability p above which a cluster is formed which 
occupies a finite fraction of the system. Let us first consider an L x L square lattice. In this 



case it is well-known that p c = 1/2 for bond percolation [9]. Alternatively, we may consider 
the ratio between the largest cluster size, sj, and the second largest one, S2 [12|. This 
measure is based on the fact that above the transition threshold, the second largest cluster 
becomes negligible with respect to the largest one, and accordingly, as soon as p exceeds 
a critical value from below, S2/S1 vanishes in the large lattice limit. Since it implies that 
one large unique cluster dominates the whole system, we call this threshold p u . Figure Ufa) 
illustrates that the ensemble average of Sifs\ scales quite well with sizes for the square 
lattice. This finite-size scaling supports the well-known result that p u = p c = 1/2 and the 
critical index v = 4/3. We will use this type of finite-size scaling methods throughout the 
paper. 

An alternative manifestation of percolation is the number b of boundary points connected 
to the middle of the lattice via occupied bonds. We will also use this concept of midpoint 
percolation throughout the paper. Around the threshold of the midpoint percolation, which 
we will call p m , a cluster penetrates the whole system and one expects a finite-size scaling 
of the form 

b = L K f[(p-p m )L^] (1) 

with k = u/(l + u) = 4/7 « 0.57 0, Q. Again, Fig. DJb) shows an excellent scaling 
collapse with k = 4/7 and p m = p c = 1/2. Consequently, the finite size scaling of b gives a 
practical alternative way of investigating percolation properties. Yet another measure, based 
on the same midpoint percolation concept, is the fraction of the boundary points connected 
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FIG. 1: (Color online) Scaling plots for regular square lattices, using (a) the ratio between the 
largest and the second largest cluster sizes, and (b) the number of boundary points connected 
to the midpoint of the lattice. The two-dimensional percolation scaling exponent v = 4/3 and 
K = + u) = 4/7 are used for scaling collapse with p u = p c = 1/2. The average is taken over 
10 6 independent realizations for each plot. 

TABLE I: Percolation thresholds and quantities used for their numerical detections. Here Si rep- 
resents the ith largest cluster size, b is the number of boundary points connected to the middle of 
the lattice, and B is the total number of boundary points. 



Threshold 


Meaning 


Observable 


Pu 


The unbounded cluster becomes unique. 


S2/S1 


Pm 


The boundary is connected to the middle. 


b 


Pb 


A finite fraction of the boundary is connected to the middle. 


b/B 



to the middle, b/AL, which becomes finite when p gets above pb- This measure will be also 
frequently used. For the planar lattices all these thresholds coincide so that we have only 
one critical threshold p c = Pu = Pm = Pb = 1/2- (The observable quantities used in the 
present paper are listed in Table [H) The crucial point in the present context is that this 
equality does not hold for hyperbolic lattices. 



III. PERCOLATION IN HYPERBOLIC LATTICES 



Suppose that n regular m-gons meet at every vertex, which we represent by a Schlafli 
symbol {m, n} 15j |. The resulting lattices are flat if we take {m, n} = {3, 6}, {4, 4}, |6, 3}, 
which together constitute the two-dimensional (2D) percolation universality class [9]. If 
(m — 2)(n — 2) > 4, on the other hand, the resulting lattice is known to have a constant 
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FIG. 2: (Color online) Poincare disk representation of the hyperbolic lattice {m, n} = {3, 7} [seven 
(n = 7) triangles (m = 3) meet at every vertex] with the total number of layers I = 6. Solid lines 
show occupied bonds in a realization with the occupation probability of p = 0.2. 

negative Gaussian curvature [3] , making a hyperbolic surface. Figure [2] shows a hyperbolic 
lattice represented as {3, 7}, where seven regular triangles meet at each vertex. Note that all 
of the triangles in this figure are congruent with respect to the metric used in this projection 



on the Poincare disk 



171 ] . As shown in Fig. [2] we construct the lattice in a concentric 



way so that the origin of the disk becomes the zeroth layer and its seven nearest neighbors 
constitute the first layer. Likewise, the second layer surrounding the first layer is composed 
of 21 vertices, and so on up to the Ith layer [I = 6 in Fig. [2]. The system is said to have a 
level I, and its size is given by N(l) = 1 + 7^. =1 [c+ — cL]/^/E with c± = (3 ± \/5)/2. For 
example, N(l = 2) = 1 + 7 (c+ ^ c - } + 7 {c +^ l) = 1 + 7 + 21 = 29. This formula shows that the 
number of lattice points increases exponentially with a distance from the origin O, yielding 
a nonvanishing surface- volume ratio in the limit iV — > oo. In other words, the total number 
of boundary points B is proportional to the system size N. The sizes of other lattices with 
different Schlafli symbols are listed in Table [III 

Since a <i-dimensional hypercube with a volume v has the surface-volume ratio oc v 
a hyperbolic lattice is usually called infinite-dimensional. In graph theory, an object with 
a nonvanishing surface- volume ratio is called nonamenable The exponential increase 

of lattice points as a function of the distance from the middle also constitute a practical 
challenge since the range of possible sizes is limited when using the finite-size scaling method. 
For this reason, we present numeric values only up to the second digit in this work. 
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TABLE II: System sizes for various structures, {m, n}, as a function of level I. Each lattice is 
grown up from a single midpoint in the zeroth layer. 



{to, n} 


N(l) 


{3,7} 






fi) 3 \ 


{4,5} 


i + ^E;=iK 2 + V3) j 


"(2 


- Vzy] 


{5,5} 




-(2= 




{6,4} 


1 + 2V2£} =1 [(3 + 2V2)' 


"(3 


- 2v^) J '] 


{7,3} 




"( 3 - 




{oo,3} 




-1 





A. Bethe lattice 



The binary Bethe lattice is the infinite binary tree where all the lattice points are equiva- 
lent [l8|. This means that the Bethe lattice lacks boundary points. Consequently, it belongs 
to the amenable class. This is in contrast to the Cayley tree (discussed in the following 
section) which includes the boundary even in the large lattice limit and hence is an example 
of a nonamenable graph. The percolation for the Bethe lattice is exactly solvable and the 
solution corresponds to the standard mean-field theory ^ . This standard mean- field theory 
describes the percolation transition for d- dimensional Euclidean lattices provided d > 6 |9j. 
However, it has limited applicability in the context of nonamenable graphs. 

The critical threshold p c is well-known since the early percolation theory formulated in 
the gelation process [jij. The point is that percolation on a Bethe lattice can be treated 
as the Galton- Watson branching process [20|. We pick up an arbitrary point as a root, 
and the set of all the points reached from it by i bonds is called the ith generation (the 
term generation will be used in trees, instead of layer). Let us denote w as the extinction 
probability that the branching process from the root is ended at some finite generation of 
the tree. For such a process, each bond to the next generation should be either unoccupied 
with a probability 1 — p, or occupied but eventually terminated with the probability pw. 
Since each vertex has two bonds to the next generation, the sum of those probabilities has 
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to be squared and w satisfies a self-consistency equation, w = (1 — p + pw) 2 , yielding 



When p > 1/2, the extinction probability is less than unity. From the percolation viewpoint, 
every vertex has one successor on average at p = 1/2, and accordingly, the cluster from the 
root vertex can be extended indefinitely. Consequently, the bond percolation of the Bethe 
lattice has p c = 1/2, at which unbounded clusters may be formed. For a general Bethe 
lattice denoted as {oo,n}, this generalizes into p c — l/(n — 1). 

An important point in the present context is that amenable graphs only have one perco- 
lation threshold. Thus for Bethe lattice the threshold p m [where a cluster from the midpoint 
(which in this case is any point because all points are equivalent) reaches a point arbitrarily 
far away] and p c [the threshold where the unbounded clusters contain a finite fraction of 
the whole system] are equal, i.e., p m = p c due to the equivalence between points. Also the 
three planar lattices ({4,4}, {3,6}, {6,3}) are amenable and consequently have only one 
threshold (compare with Sec. [IT])- As will be discussed in the following, nonamenable graphs 
such as hyperbolic lattices have two distinct percolation thresholds. 

B. Cayley tree, {oo,3} 

The Cayley tree is a tree grown from the root vertex up to the Ith generation where the 
root vertex is identified as the middle of the lattice. This is an example of a nonamenable 
graph. While it has sometimes been presumed that the Bethe lattice is an adequate limiting 
case of the Cayley tree, more recent studies suggest that the Cayley tree in the limit of 
/ — > oo has different critical properties from those of the Bethe lattice |l| . Note that vertices 
are not equivalent for the Cayley tree so that one can clearly define which generation a vertex 
belongs to. The branching-process argument is again applicable and the cluster from the 
root vertex reaches the bottom of the tree at p m = 1/2. On the other hand, the uniqueness 
threshold is located at p u — 1 [6j, since s 2 /si remains finite at any p < 1 [Fig. E(a)]. The 
vanishing of S2/S1 as p u = 1 is approached can be obtained as follows: The total number of 
bonds in the tree with a level I is K — 2 l+1 — 2. Suppose that precisely one bond is broken 
on average. The occupation probability corresponding to this is p = (K — 1)/K- Suppose 



w = 
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further that the broken bond connects the zth and (i + l)th generations. The probability to 
select this bond is given by 

O'i+l 

P(i) 



2'+* -2' 

Breaking one bond creates precisely two clusters. The smaller one is below the broken bond 
and has a size of s 2 = 2 — 1. The size of the larger one is consequently s\ = N — s 2 = 
2 l+1 — 2 l ~ l . The expectation value of the ratio S2/S1 is obtained as follows: 

1-1 

s 2 \ s 2 



si f si 



2 l ~* - 1 2 i+l 
1^ 2 l+1 - 2 l ~ l 2 l+l - 2 



i=0 

1-1 



2 1 - 1 - 1 2 i+1 
2-^< 2'+i _ 2 l ~' L 2 l+1 — 

i=0 

1 l ~ l 1 

— Y - 



i=0 



Since p = (if — 1)/K, one can express this result directly in terms of p. By using the 
connection 1 — p — K^ 1 pa 2~ l ~ 1 , one obtains 

" x oc -(1 -p) log(l -p). 



,81, 

Thus the ratio S2/S1 vanishes as —(1 — p) log(l — p) as the threshold p u = 1 is approached 
from below. This is illustrated in Fig.[3](a): The Cayley tree has two thresholds and p m < p u . 
The midpoint percolation quantity b instead possesses the scaling form 

b = rf[(p- Pm )i 1 /% (2) 

with k = and v = 1. This scaling form can be derived as follows: The number of boundary 
points which is reached from the midpoint for a given value p is b = 3p x (2p)' _1 from which 
it follows that 

3 3 
b = - exp[log(2p)Z] pa - exp[2/(p - p m )\. (3) 

The validity of this finite-size scaling is demonstrated in [Fig. EJ^b)] together with the exact 
scaling form given by Eq. ©. 
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FIG. 3: (Color online) Numerical results in Cayley trees, averaged over 10 6 realizations, (a) 
The ratio between S2 and s\. Inset: A closer look near p = 1, where one can hardly see any size 
dependency, (b) Scaling plot by Eq. ((2J). The scaling function is well-described by | exp[2/(p— p m )\ 
near p m = 0.50. (c) as a function of p, provided that b(p) ~ B^. The horizontal axis is in the 
log scale, confirming = 1 + log 2 p- (d) Cluster size distribution at various p with I = 15. It 
approaches to P(s) ~ s~ 2 as p —>■ 1, while finiteness of the system adds an exponential cutoff. 

One may also note that since I ~ log 2 N, it follows that b ~ (2p) z = Np l pa Np log2 N = 
A rl+log2 p . This means that b scales as iV^ (or equivalently as B^) at any value of < p < 1 
with a p-dependent exponent 0=1 + log 2 p [Fig. 12(c)] . A direct consequence of this is 
that b/B eventually goes to zero at every p < 1 in the limit of infinite N. Thus b/B 
is discontinuous at pb — 1 for / = oo. This is in contrast to the ratio s^/sx which goes 
continuously to zero at p u = 1. 

One can also obtain the limiting behavior of the cluster size distribution P(s) by removing 
one bond. In this case P(s) ds is the probability of finding a cluster with a size in the interval 
[s, s + ds] and P(i) di is the probability of the single bond to be broken between generations 
i and (i + 1). This latter probability is also equal to the probability of finding a second 
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largest cluster of size s 2 oc s L \ Thus P{s%)dsi = P(i)di from which follows that 

di 



P(s 2 ) = Pit] 



ds-: 



oc s 2 2 , 



where the last proportionality is obtained from s 2 oc s L_ \ This suggests that the size 
distribution of clusters P(s) (excluding the largest one) should approach the form P(s) oc s~~ T 
where r — > 2 as p — > 1. Since the generation i in a tree is actually a discrete variable it 
follows that the possible sizes of S2 are also discrete and this discreteness becomes noticeable 
as p — > 1. This is illustrated in Fig. EJ^d) which shows the approach to P(s) oc s~ 2 as well 
as the log-periodic oscillations caused by the discreteness. 

A more intuitive explanation for P(s) may be gained by mapping the percolation problem 
in a Cayley tree to a branching process, namely, the formation of family trees [2lj |: Let us 
consider the family-size distribution in the z-ary tree (z = n — 1) at the fcth generation, 
with the net birth rate A = log(zp) as each family grows by (zp) k . When a bond is broken 
in a tree graph, the top of this detached branch is interpreted as the first ancestor of a new 
family, and we describe this top vertex as an immigrant, whose number is proportional to 
the existing population size A^. If the number of immigrants is written as (N^, the birth 
and immigration should yield a constant growth of population, £ + A = log z, as we know 
Nk+i = zNk- According to Ref. [211 ]. the family size distribution at the fc-th generation 
exhibits a power law with an exponent r' = 2 + Q/\ = 2 + log(l/p)/ log(zp) as k — > 00. 
Since the volume of a given tree is proportional to its surface, it seems plausible that r ~ r' 
for the cluster size distribution, P(s). Note that r' = 2 for p = 1 which is also the limiting 
result for the Cayley tree. This gives a hand-waving argument suggesting that the size 
distribution form P(s) oc s~ T could be valid also for some range of p below one. 



C. Heptagonal lattice, {7,3} 

Now we consider a heptagonal lattice, denoted as {m, n} = {7,3}. Since the probability 
to find a loop is roughly 0(p m ), we expect the results obtained for a Cayley tree will remain 
valid to some extent. The lower threshold is determined from the finite-size scaling of b 
and Eq. (j2J) as shown in Fig. HJ^a). This determines the value of p m ~ 0.53. We note that 
this is rather close to the exact result for the Cayley tree l/(n— 1) = 1/2, and shows the 
same scaling form as in Eq. (j2j). This suggests that the loops have only a small effect on 
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FIG. 4: (Color online) Numerical results for heptagonal lattices, {7,3}, averaged over 10 6 realiza- 
tions, (a) Measurement of b yields a consistent result with the tree case, (b) The ratio between two 
largest cluster sizes S2/S1 gives p u « 0.72, and (c) extrapolating b/B gives pb ~ 0.72, supporting 
Pb = Pu- (d) The exponent 4> as a function of p shows a clear difference from the tree case [see 
Fig. [3|c), for comparison]. At p ^ 0.84, <j) largely fluctuates as it is hard to determine the size 
dependency from the numerical data. 

the percolation properties at small values of p, as is also suggested by the actual structure 
of the clusters illustrated in Fig. EJ 

The upper threshold is determined from the finite-size scaling of S2/S1 which gives p u ~ 



0.72. According to a dual-lattice argument 



22] , a percolation threshold for a regular lattice 



{m,n} is predicted approximately asm/(m + n). Applying this argument somewhat ad hoc 
to the present case yields 7/(7 + 3) = 0.7, which is in fact fairly close to the actual value, 
Pu « 0.72. 

Alternatively, the upper threshold can be determined from b/B by assuming that the size 
scaling form is the same as for the Cayley tree, 

b/B ~ dA^- 1 + c 2 , 

and extrapolating the numerical results from / = 6, . . . , 11 to the infinite-size limit [Fig.[4](c)]. 
This extrapolation gives \b/B\ < O(10~ 3 ) below p < 0.72 and becomes positive finite above 
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FIG. 5: (Color online) Scaling plot of b/B in {7,3} using Eq. (@J, with p b = 0.72, = 0.82, and 
l/z7 = 0.12. 
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FIG. 6: (Color online) (a) Cluster size distribution, P(s), at various p in {7,3}. It shows a 
power-law form also above p u . (b) A schematic view of the Poincare disk centered at O. The 
lattice is constructed up to the dashed circle, which is apart from the circumference of the disk by 
AB = CD = 5 in the Euclidean measure. 

that, suggesting pb ~ 0.72. It thus suggests pb = p u within our numerical accuracy. This 
corresponding scaling exponent <j)(p) is plotted in Fig. IU(d). While <ft(p) is an increasing 
function of p in the Cayley tree, it is a convex function in {7,3}. The crucial difference is 
that is still less than one at p = pb- In order to study the transition at pb, we use the 
finite-size scaling form: 

b/B-N^fiip-pb)^*}. (4) 

The scaling collapse at pb = 0.72 determines the critical indices to <fi ~ 0.82 and v « 0.12 
(Fig. W>. 

Another crucial difference from the Cayley tree is the appearance of a supercritical region 
at p > p u . We find that the hyperbolic lattices in this region display a power-law behavior 
in the cluster size distribution [Fig. [6]^a)]. A hand- waving argument for this behavior goes 
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FIG. 7: (Color online) Numerical results for {4,5}, averaged over 10 6 realizations, (a) Scaling plot 
with p m = 0.27. (b) Scaling collapse with = 0.52, (ft = 0.82, and \jv = 0.11. 

as follows: The probability of finding a cluster with a certain size is dominated by its surface 
area, since connections through the surface have to be cut to isolate this cluster from the 
surrounding. It is thus believed that -P(s) oc exp[— rj{p)s 1 ~ 1 l d ] in a d- dimensional lattice with 
some p-dependent constant r](p) j^l, Q. For a hyperbolic lattice with d = oo, one may well 
expect P{s) oc e~ s . In case of our nonamenable setting, however, most large-sized clusters 
are facing the outmost boundary of the lattice, where the outward connections are already 
absent. Suppose a cluster contained in a fan-shape OBC on the Poincare disk, where the 
lattice is constructed up to the dashed line in Fig. [6](b). If AB = CD = 5 in the Euclidean 
measure, the hyperbolic length of the arc BC connecting B and C is of an order of <5 _1 . 
Since the cluster size is closely related to the surface length in a hyperbolic lattice, which 
is proportional to [17], we may say that s oc 5~ l . At the same time, the hyperbolic 
length that one should cut out to isolate this cluster is roughly the length of the geodesic 
connecting B and C, which grows as log <5 _1 . Note that BC needs not be considered here 
since it is a part of the lattice boundary, and this makes the fundamental difference from 
the exponential decay. In other words, the number of bonds cut for hyperbolic lattices is 
not proportional to the size s, but only to log s. This gives the cluster size distribution 
in a power-law form as P(s) ~ exp[— r/(p) log s] = s~ v ^ p \ It is also possible to infer that 
7](p) should be an increasing function of p, as clusters are merged to the largest one in the 
supercritical phase. 
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FIG. 8: (Color online) Numerical results for {3, 7}, averaged over 10 6 realizations. Here are shown 
the scaling plots (a) at p m = 0.20 and (b) at pb = 0.37 with <p = 0.82 and 1/v = 0.11. 

D. Comparison with {4, 5} and {3, 7} 

In order to investigate the generality of our results, we also study two additional hy- 
perbolic lattices with different structures {m,n}. In all cases we find precisely the same 
scenario. For the hyperbolic lattice {4, 5}, we find p m ~ 0.27, p u = pb ~ 0.52 [Figs. [7(a) and 
UJb)]. The scaling behavior at the lower threshold p m is the same as for the Cayley tree and 
for the lattice {7,3}. At the upper threshold pb = 0.52, the critical indices <fi = 0.82 and 
1/v = 0.11 was determined from the scaling collapse [Fig. Etc)] . These values are identical 
to the ones found for lattice {7, 3} within numerical accuracy. 

For the lattice {3, 7} which is dual to {7, 3} the same agreement is found: The same size 
scaling is found at the lower threshold p m ~ 0.20 [Fig. Eta)]. At the upper one, p u ~ 0.37, 
the critical indices = 0.82 and 1/v = 0.11 are found [Fig. E(b)]. This is again in striking 
agreement with the lattices {7, 3} and {4, 5}. Thus our results are consistent with a universal 
critical behavior at the second threshold for all hyperbolic lattices {m, n} provided both m 
and n are finite numbers. This critical behavior is distinct from the tree case {oo, n}, which 
has = 1 at p = pb- The present accuracy suggests that the critical indices are to good 
approximation »s 0.82 and 1/v « 0.11 ± 0.01. 

We also note that while the lower threshold p m is still close to the tree result l/(n — 1) = 
1/4 for {4,5}, the deviation becomes large for {3,7}, and that the estimate Pb = Pu = 
m/(m + n) in general only gives a very crude estimate of the upper threshold. 



14 



IV. SUMMARY 



We have investigated the percolation thresholds and the critical properties of the per- 
colation transitions for hyperbolic lattices using finite-size scaling methods. Two distinct 
percolation thresholds were found: The lower one corresponds to the threshold when the 
probability of finding a cluster from the midpoint to the boundary becomes finite and the 
second when the cluster containing the midpoint, with finite probability, also contains a 
finite fraction of the boundary. This is in contrast to the planar lattices which only possess 
a single percolation threshold because the two thresholds above coincide. The Cayley tree 
was used as a benchmark. It was found that the lower threshold for the hyperbolic lattices 
has the same scaling properties as for the Cayley tree and that the power-law dependencies 
characterizing the region between the two thresholds are also like Cayley trees. However, 
the second higher threshold has a different critical behavior. Our results are consistent with 
a universal behavior at the higher threshold for all hyperbolic lattices {m, n} with m and n 
finite. This critical behavior is characterized by two critical indices 0.82 and v 0.11. 
What actually determines these critical indices is still an open question and will be the 
subject of future research. 
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